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*0 ' Abstract 

cn : 

| We study the impact of higher dimension operators in the inflaton Lagrangian on the non-gaussianity 
■ of the scalar spectrum. These terms can strongly enhance the effect without spoiling slow-roll, though 
it is difficult to exceed /nl ~ 1, because the scale which suppresses the operators cannot be too low, 
—, , if we want the effective field theory description to make sense. In particular we explicitly calculate 
the 3-point function given by an higher derivative interaction of the form (V0) 4 , which is expected 
| to give the most important contribution. The angular dependence of the result turns out to be quite 
different from the minimal case without higher dimension operators. 
^ ' 

P ■ 1 Introduction 

£3 ■ 

In most models, inflation takes place at energies many orders of magnitudes higher than those 
accessible with accelerators. It is therefore mandatory to explore all what we can learn about 
^ . this high energy regime from the signatures left by inflation in the present Universe. In the 
effective field theory framework the effects of short distance physics can always be encoded 
in higher dimension operators for the inflaton Lagrangian 1 0. The obvious problem is how 
to experimentally probe the presence of these operators. Despite the successes of inflationary 
cosmology in fact, there are very few observables which can be used to constrain the inflaton 
dynamics. 

References [2] and [3] studied the possibility that higher dimension operators can change 
the prediction for the tilt of the gravitational wave spectrum. This quantity is fixed, in single 
field models, once the ratio between scalar and tensor modes is known, through the so-called 
'consistency relation'. The problem with this kind of observable is that the detection of 
gravitational waves is possible only if the scale of inflation is close to the presently maximum 
allowed value, while most models predict a totally negligible production of tensor modes. 
Even if we restrict to models with a sensible production of gravity waves, the tilt of this 



1 Another logical possibility is that unknown high energy physics modifies the vacuum state of the inflaton 
in the approximately de-Sitter background of inflation pp. We will not consider this possibility here, assuming 
that the inflaton is in the usual Bunch-Davies vacuum. 
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spectrum remains almost impossible to probe. Order one deviation from the consistency 
relation could be measurable with dedicated polarization experiments as CMBpol, only if the 
tensor to scalar ratio is very big: T/S > 0.1 |3]. 

The purpose of this note is to indicate a potentially better test of higher dimension op- 
erators in the inflaton Lagrangian, which is not related to the elusive tensor modes: non- 
gaussianities in the scalar perturbations. The basic reason why this could be a good smoking 
gun of high energy effects in inflation is that it is quite suppressed in the standard framework. 
This follows from the fact that the inflaton must be very weakly coupled to have a slow-roll 
phase, so that its fluctuations are very close to those of a free field, i.e. gaussian. 

Consider for example the effect of the cubic self-interaction of the inflaton 0, given by 
V^ 3 \(f)), where V is the inflaton potential. The amount of non-gaussianity is parametrized in 
this case by V^((p)/(2TiH), with H the Hubble constant during inflation. In the equation of 
motion for the inflaton fluctuations, this ratio gives the relative importance of the non-linear 
term with respect to the free field terms at horizon crossing. The slow-roll conditions impose 
a strong bound on the derivatives of the potential; in particular the constraint on the third 
derivative is 

M> p — «1 _«___ = 3 ^ with ^jM|(- 

The limit is quite strong and it is completely model independent, because P-n is the measured 
curvature power spectrum: 3P^/ 2 ~ 1.4 • 10~ 4 . We therefore see that the existence of a 
slow-roll phase requires that the inflaton fluctuations are very close to gaussian. 

The full calculation of the 3-point function has been carried out quite recently in [5] and 
jHj- It turns out that the strong constraints we described on the inflaton potential make the 
inflaton coupling to gravity the leading source of non-gaussianity. The result is suppressed 
with respect to (Q) by a combination of the slow-roll parameters. Parametrically the level of 
non-gaussianity is thus predicted to be 

0(e, V ). J-^LA< 3 .10- 6 W ith V = M P ^. (2) 

Z7T y/e Mp V 

Comparing this prediction with the current limit coming from the WMAP satellite [7] , which 
is of order 5 • 10~ 3 , we see that the level of non-gaussianity in conventional inflationary models 
is quite low and certainly it will not be measured in near-future experiments. This very low 
prediction is clearly good to test deviations from the simplest models. 

A possibility to get sensible non-gaussianities in to have an additional light field {i.e. with 
mass smaller than H) during inflation. This new degree of freedom is not constrained from the 
slow-roll requirements and if its fluctuations are somehow converted into metric perturbations, 
these can be much less gaussian than in the usual scenario. This possibility is realized in 
multiple-field models of inflation [HJ E], in the curvaton scenario ^U] and in the recently 
proposed mechanism with curvature perturbations created by fluctuations of the decay rate 
of the inflaton [TT1IT2]. 

We will follow a more minimal possibility here. We want to see if the addition of higher- 
dimension operators in otherwise conventional models can give an increase of the result ((21). 
Before concentrating on an explicit calculation in section El we discuss which operators we 
are allowed to add in an effective field theory description of inflation. 
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2 Effective field theory during inflation 



During inflation there are two energy scales of interest. One is related to the classical time 
evolution of the inflaton and it is given by 1//2 . The other one is typical of the quantum 
fluctuations of the inflaton around this classical background and it is given by the Hubble 
scale H. The separation of these two energy scales is fixed by the normalization of the scalar 
perturbations originally measured by COBE: if/0 1 / 2 ~ 10 -2 . As we want to encode unknown 
short distance physics in an effective field theory description, it is natural to ask which is the 
typical scale of this theory or equivalently which heavy states we are allowed to integrate out. 

In principle, as we are interested in the study of the inflaton perturbations, which are 
characterized by H, we could be tempted to write down a generic Lagrangian for the fluc- 
tuations, with higher dimension operators suppressed by a scale which is only required to 
be bigger than H. This is for example the attitude taken in This procedure is too ar- 
bitrary. Many constraints on the Lagrangian come from the requirement of having a slowly 
rolling background or, in other words, it is not easy to imagine that any effective operator 
for the inflaton perturbations can be obtained by a UV completion compatible with slow-roll. 
Moreover if we took this point of view, there would be no reason to expect nearly gaussian 
perturbations at all, because we could write down a generic self-interaction for the inflaton, 
while we saw that the constraints come from the required flatness of the potential in order to 
have an inflating phase. 

Our approach will be to consider a generic effective Lagrangian, which encodes all the 
heavy physics above the higher scale 1 / 2 . In principle we are allowed to write down all 
operators we want, with the only constraint that an inflating phase is possible. Every choice 
will modify the background evolution and the fluctuations at the same time. 

The question now is which operators can give big non-gaussianities without spoiling the 
slow-roll of the inflaton. References jSJ E] conclude that for a generic potential we get the 
estimate @. We are therefore led to consider terms with derivatives. It is not difficult to 
conclude that operators with derivatives acting on each field 0, i.e. compatible with a shift 
symmetry — > + c, are the best candidates because they do not interfere with the slow-roll 
evolution of the inflaton. In particular one could impose an approximate shift symmetry on 
the inflaton, so that only these terms are allowed in the limit of exact symmetry. Obviously 
the symmetry cannot be exact otherwise the inflaton would be a Goldstone boson with a 
perfectly flat potential; some tilt must be introduced through an explicit breaking, making 
the inflaton a pseudo Goldstone boson 2 . 

It is therefore interesting to study the set of all operators compatible with the exact shift 
symmetry of the inflaton 3 , suppressed by the appropriate power of the scale of new physics M. 
These derivative interactions could arise, for example, from a Higgs UV completion, where 
the shift symmetry is derived from a linearly realized, spontaneously broken U(l) invariance, 
or they could be the consequence of strong dynamics at the scale M. We stress that these 
operators, modifying just the kinetic part of the inflaton equation of motion are not dangerous 
for the slow- roll and we are allowed to take the scale M to be as low as 1//2 : when this limit 
is reached all the higher-dimension operators are equally important for the classical evolution 

2 The idea of the inflaton as a pseudo Goldstone boson was introduced a long time ago in ^3] and recently 
reinterpreted in the context of higher-dimensional (or higher-dimcnsionally inspired) ideas |14|. 

3 The effects of higher derivative interactions on the consistency relation have been studied in EJ. 
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of the inflaton and the effective field theory description ceases to be useful 4 . 

Among all the operators compatible with the shift symmetry, the most important are 
those of the form (V0) 2 ™. Operators with more derivatives acting on one field (as □</>• (V0) 2 ) 
are less relevant because each derivative gives a suppression of order H/M, both in studying 
the evolution of the classical background and for the fluctuations around it. In particular the 
lowest dimension operator will be (V0) 4 /M 4 . 

What is the size of the 3-point function we expect adding this higher- dimension operator? 
One of the four legs will be evaluated on the classical background ~ V'/H, so that we 
obtain an interaction among three fields with derivatives acting on them; each derivative 
gives a factor of H at horizon crossing. The corrections to the free theory, which control the 
deviation from exact gaussianity, are thus given by 

V H 3 1^1 V'H 1 1 H 

— ■ ■ ~ ( 3 1 

H M 4 2nH ~ 2tt eV 2vr M P ' V; 

where the upper limit comes from setting M 4 = 2 ~ e V ( 5 ). As we discussed, in this limit 
all the higher terms become important and the effective field theory description breaks down. 

We see that the effect can be much bigger than the result obtained in the minimal case 
(J2J), which is further suppressed by a slow-roll parameter. The amount of non-gaussianity is 
still small (< 5 • 10 -5 ) and, as we will discuss, it is close to the limit of what is measurable. To 
make more precise statements about the size of the effect and to study the angular dependence 
of the 3-point function we now turn to the explicit calculation. 



3 Calculation of the 3-point function 



The action we will consider is 



S 



d xa 



-M 2 P R 



-(V0) 2 -n0) + ^r7i(V0) 2 (V0) 2 ... 



(4) 



2~-r~ 2 v ■ r, ■ vr, ■ gM 4 

where the dots stand for additional higher dimension operators compatible with the shift 
symmetry — > + c. The background solution for the metric is of the form 

ds 2 = -dt 2 + e 2Ht dx i dx l , (5) 

with H a slow- varying function of time (H/H 2 <C 1). From the action (pjj it is easy to get 
the equations describing the evolution of the scalar and H [T^l E] 



H 



1 



3 0^ 
2M 4 



W) 2 r 8M 4 



+ 3#0 



1 + 



2M 4 



+ V'{4>) = . 



(6) 



(7) 



4 We stress that it would be possible to add specific operators with many derivatives (as D(f> ■ (Vcj)) 2 for 
example) suppressed by a lower mass scale to get big non-gaussianities without spoiling slow-roll. On the 
other hand, without an explicit UV completion, we are forced to take all operators compatible with the given 
symmetries, suppressed by the same scale of new physics M and to remain in the regime in which this effective 
field theory description makes sense. It is conceivable, but unlikely, that a specific UV completion gives a 
much bigger non-gaussianity than what obtained here. 

5 In the same limit we obtain order one deviations from the consistency relation 
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As we are interested in the calculation of the 3-point function at first order in 2 /M 4 , we will 
often neglect the 1/M 4 corrections to eqs. © and (J7J). In this limit the background equations 
reduce to the usual ones describing single-field inflation. 

To study the perturbations around the background solution we follow f5|, using the ADM 
metric 

ds 2 = -N 2 dt 2 + h ij (dx i + N i dt)(dx j + N j dt) . (8) 

In this formalism the dynamical degrees of freedom are and the scalar 0, while N and N l 
have no time derivatives in the action, so that they behave as Lagrange multipliers. 

We now have to fix the gauge to second order in the fluctuations. As discussed in jS], 
there are two useful choices. The first one is given by 6 

5<p = , h i:i = e 2Ht+2 %j , det h = . (9) 

The variable C is the non-linear generalization of the one introduced by Bardeen, Steinhardt 
and Turner [T^]. In references ^7] and [5], it is proven that C is constant outside the horizon. 
For wavelengths much bigger than the Hubble radius, we can neglect all gradient terms, so 
that different parts of the Universe follow the same evolution, with a certain shift in time 
parametrized by (; this explains intuitively why this variable stays constant outside the 
horizon. Note that the same result holds true in our case, after the addition of the higher 
derivative terms: these interactions always contain derivatives acting on the fluctuations 
around the background, so that they are irrelevant outside the horizon. 
The other useful choice of gauge is given by 

= (p{t) + <p(t, x) , = e 2H %j , det h = . (10) 

This gauge has an important property that will be useful for our calculation. At leading order 
in the slow-roll parameters the action for tp can be calculated from (0J), without considering 
the gravitational contributions contained in y/—gR [5J. This is quite intuitive, because in the 
limit of a completely flat potential (p ^ does not give any spacetime curvature. 

In this gauge the three point function (<p<p<p) can therefore be calculated at leading order 
in the slow-roll parameters from the interaction Lagrangian 

3Ht 

Ant = -|^(#)(-0 2 + e- 2H \d lV f ) , (11) 

where one of the four legs of the (V0) 4 interaction in (jlj) is evaluated on the time-dependent 
background <fi. As discussed we can disregard, at leading order, the 1/M 4 corrections to the 
background in eqs. © and ((7j). 

Before calculating ((pcpcp), we have to discuss an important point. We are interested in 
the (CCC) correlator: this remains constant outside the horizon and reappears as a non- 
gaussian contribution to the perturbations we observe today. Once we cross the horizon we 
have therefore to go from the ip variable to ( and this must be done at second order in the 
pertubations: a second order contribution could give an additional piece (CCC) ^ (w) iff) ■ 

6 As we are not interested in the tensor modes, we will not specify the parametrization of hij. See 
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This is what happens in the calculation [S] without higher derivative terms; the relation 
between ( and ip is schematically given by 

( = -jV + 0(e, V )(j V y . (12) 

The non-linear term is proportional to the slow-roll parameters and it exactly cancels the p 
evolution outside the horizon to give a time- independent (. We will neglect this non-linear 
term in our calculation as it is sub-leading in the slow-roll expansion. Moreover the relation 
((p>) is not modified outside the horizon by the addition of the higher dimension operator 
(V0) 4 /M 4 : as we discussed, this term contains derivatives acting on the fluctuations and it is 
therefore exponentially irrelevant when the wavelength exceeds the Hubble radius. So we can 
take just the linear term in (fT2*|) to convert the result in the ( variables after horizon crossing. 

We now turn to the explicit calculation 7 of (tpp>p>). At leading order in slow-roll we can 
take H constant and we can write the contribution to (ipipip) from the interaction (jllj) as an 
integral over the conformal time r\: 

fat. wpt> = ^ ( E *<) n<sf) IS- < 13 > 

f° - - 

■i dr](- k\k\k\ rf - (k x ■ k 2 )kl(l - ifar])(l - ik 2 r]))e lktV + perm. + c.c. 

J — oo 

where fcj are the moduli of the wavevectors and k t = k\ + k 2 + & 3 . The integration is done 
with an infinitesimal contribution in imaginary time (r] — > 7]+ie\r]\); this prescription projects 
the initial state in the usual Bunch-Davies vacuum at past infinity. 

Note that the perturbations (p are strongly coupled well inside the horizon when their 
typical wavelength is of order M: all the higher dimension operators become important. 
Does this imply that our results depend on the UV completion of the theory? Obviously no. 
Whatever the UV completion might be, the evolution of the background is very slow with 
respect with the typical scale of this unknown physics because M ^> H: we therefore expect 
that, following the slow expansion, the system is always in its vacuum state. Although we do 
not know how this state looks like for very short wavelengths, when the description in terms 
of a scalar (p breaks down, we know that below the scale M the system is in the Bunch-Davies 
vacuum of the weakly interacting (p field. The same situation is at the origin of the so called 
'trans-Planckian problem': the description of the inflaton as a weakly interacting field ceases 
to make sense at the Planck scale. Even if the details of quantum gravity are unknown, it is 
likely that all the predictions are independent of this regime, as in our simple example. 

Returning to the calculation, after the transformation to the ( variable through £ = — S^p 

we get the final result performing the integral in (fTH|) 

< i4 > 



7 For all the details about this kind of calculation we refer the reader to the very clear discussion of 
whose notation we have mostly followed. 
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As one could have expected the angular dependence given by the second line of (fl"4"j) is 
completely different both from the result of the minimal case jH] and from that assumed in the 
analysis of the data [7] which follows from the assumption that the non-gaussianities come 
from a field redefinition 

c = c g - ||/ni4 2 , (15) 

with ( g gaussian 8 . The most relevant difference is the behavior in the limit in which one of 
the wavevectors goes to zero. In this limit one of the fluctuations has a very big wavelength, 
it exits the horizon and freezes much before the other two and acts as a sort of background. 
What do we expect in this limit? Let us take k 3 very small and consider the case in which 
a spatial derivative diip in the interaction Lagrangian (JTTjl is acting on (p%. The 2-point 
function (pi(f2 depends on the position on the background wave and it is proportional to diip 3 
at linear order. This variation of the 2-point function along the ip 3 wave is averaged to zero in 
calculating the 3-point function { 1 4>% X ( P% 2 ( P%^) , because the spatial average (ip 3 diip 3 ) vanishes. 
So we are forced to go to the second order and we therefore expect that the second line of 
(|14|) behaves as k\ and not as ^3 in the limit k% — > 0. Note that this result holds also when a 
time derivative acts in (jllj) on the long wavelength mode, because the scalar wavefunction in 
de-Sitter space goes to a constant outside the horizon as the square of the physical wavevector 
fc|/(aif) 2 . Although it is not evident, it is easy to check that the second line of ([Hj) does 
indeed go to zero as k\ when fc 3 — > 0. This is a radical difference from the behaviour studied 
in [5] and [7j: they both go to a constant in this limit. The difference is a consequence of 
the derivative interaction (V0) 4 , which favors the correlation among modes of comparable 
wavelength. 

The different angular behaviour will probably change the limits coming from data analysis. 
It is anyway interesting to calculate the "effective" /nl in (|15|) in the case of an equilateral 
triangle: k\ = k 2 = k 3 . It is easy to get 

,equil. _ / 16 n 

The result is not far from the naive estimate we did in the introduction and it should be com- 
pared with what is obtained in absence of higher derivative terms: in the same configuration 
ref.s |S] and [Hj obtain 

Comparing (|TB|) and (jl7|) we see that for sufficiently low mass scale M the new contribution 
dominates over the usual one, which is suppressed by the slow-roll parameters. As we discussed 
it is not possible to take M arbitrarily small, because at a certain point the expansion <p/M 2 
ceases to make sense and we should consider all the higher dimension operators in the action 
(@J). We can estimate the lower limit on M assuming that there is some strong coupling regime 
around M and using naive dimensional analysis to estimate the size of the operators. In this 
case we expect a Lagrangian given by 

A1 ^^(v^ + ^(v0) 4 + l-L(V0) 6 -- 



16tt 2 V2A 2V ^ 8A 4V ry 48 A 6 



8 This pattern of non-gaussianity which is local in real space is characteristic of models in which the non- 
linearities develop outside the horizon. This happens if the perturbations are created by fluctuations in the 
reheating efficiency an d in the curvaton scenario 
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In front of each of operator unknown order one coefficients should be understood. Going 
from the dimensionless scalar <fi to the canonical normalization we obtain that the expansion 
makes sense for (p/M 2 < 1. The same conclusion holds if we take a UV completion where 
represents the Goldstone boson of a spontaneously broken U(l) symmetry. 

Even if we could have some order one enhancement, e.g. if the (V0) 4 operator is somewhat 
bigger than its natural size, we conclude that a reasonable limit on non-gaussianities in single- 
field models is 

/nl < 1 ■ (19) 

This limit is roughly what is expected to be measurable with the forecoming CMBR 
experiments. With the new data WMAP will reach a sensitivity |/nl| ^ 20, while the Planck 
satellite will improve the limit to |/nl| ^ 5 [THj. Note that the latter result is not far from 
what can be measured in principle in a CMBR experiment. The number of independent 
spots in the sky is limited by the diffusion length of photons before recombination and Planck 
sensitivity will not be far from this limit. Large scale structure surveys are not so severely 
limited by cosmic variance, but it is seems difficult that they can reach limits on primordial 
non-gaussianities stronger than those coming from CMBR experiments ^H]- The limit ([19)1 
is therefore at the threshold of detectability. 

We have to stress that many additional effects are expected to give /nl ~ 1, even starting 
from a gaussian initial perturbation. All the non-linearities in the evolution from the im- 
printing during inflaton to the detection of temperature fluctuations will give /nl ~ 1, as it 
is clear from the definition (fT3j) . These effects should be under control before a clear signal 
of primordial non-gaussianity can be found. The characteristic angular behaviour ()14|) can 
be a powerful discriminator against all the other effects. As the 3-point function contains in 
principle much more information than the single parameter /nl, it would be very interesting 
if the data analysis were carried out considering different angular patterns [2Uj . 

4 Conclusions 

Due to the strong suppression of non-gaussianities in minimal models of inflation, a deviation 
from this prediction is probably the most efficient method to detect the effect of short distance 
physics during inflation, which can change the inflaton dynamics through higher dimension 
operators. 

We studied the effect of operators compatible with a shift symmetry of the inflaton as 
the most promising set of higher dimensional terms. The explicit calculation of the three 
point function (fUjl shows a particular angular dependence as a consequence of the derivative 
interactions. Although the size of the operators can only be calculated once the UV completion 
is known, we conclude that /nl ~ 1 represents the upper limit of this effect, before the 
effective field theory description breaks down. This limit is close to what will be measured 
by the Planck experiment. 

This result can be considered in some sense an upper limit for single field models. A 
detection of a much bigger non-gaussianity, f^L 3> 1, would point towards models with 
additional light degrees of freedom during inflation. 
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